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ABSTRACT 

Safety  analyses  of  a  loss-of --coolant  accident  in  a 
pressurized-water-reactor  nuclear  power  plant  require 
modeling  the  flow  through  the  main-coolant  pumps  in  any  one 
of  several  extreme  operating  conditions.   A  semi-empirical 
method  of  predicting  the  two-phase  performance  of  centri- 
fugal pumps  has  been  proposed  and  was  applied  to  normal  pump 
operation  with  good  results.   In  this  report  the  proposed 
method  is  applied  to  the  reverse-flow,  reverse-rotation  mode 
of  operation,  also  with  good  results. 

By  making  simplifying  assumptions,  the  flow  through 
the  pump  is  idealized  and  expressions  for  the  theoretical 
pump  characteristics  derived  from  Euler's  equation  for  the 
mean  flow  through  a  turbomachine .   General  expressions  for 
the  flow  losses  in  two-phase  flow  are  derived,  and  the  flow 
losses  expressed  as  the  difference  between  the  theoretical 
characteristics  and  the  actual  performance  of  the  pump. 
The  ratio  of  the  losses  in  two-phase  flow  to  those  in  single- 
phase  flow  is  shown  to  be  primarily  a  function  of  pump 
geometry  and  fluid  void-fraction.   Using  experimental  data 
from  single-  and  two-phase  flow  tests  on  a  pump  of  known 
geometry,  the  ratio  of  the  losses  was  calculated  for  various 
conditions  of  flow  and  plotted  as  a  function  of  void  fraction. 
The  results  plotted  in  this  fashion  show  excellent  correlation 
and  the  function  can  be  used  to  predict  the  two-phase 
characteristics  of  similar  pumps. 
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LIST  OF  SYMBOLS 


A       flow  area 

a       two-phase  flow  function  e  [-r^—] 
r  1-a 


p 
a  ,   v 


pL 


c       fluid  velocity,  absolute 

f       two-phase  flow  function  (eq.  8a) 

g       gravitational  acceleration 

q       constant  in  Newton's  Law 
^c 

H       head 

AHot  th  ~  AHot 
H*      head-loss  ratio  =   -rrz — *- 777 — ° 

An        J.-U        ~        All      _ 

ospth     osp 
h       enthalpy 

K,  K.  ,  Kj./  Kfi         constants  defined  by  eqs.  12,  13,  and  15 

• 

m  mass  flow  per  unit  time 

N  rotor  speed,  rev/min 

p  fluid  pressure 

Q  volume  flow  per  unit  time 

r  impeller  radius 


R 


1-x  2 

friction-factor  multiplier  =  N ] 

^        1-a 


s  slip  velocity  ratio  =  (C  /CT ) 

U  rotor  peripheral  velocity 

w  fluid  velocity  relative  to  rotor 

x  quality  =    (mv/mT) 

a  void  fraction  =  (A  /A  ) 

3       angle  of  fluid  vector  relative  to  tangent  to 
rotor  periphery 

6'      angle  of  tangent  to  rotor-blade  centerline  relative 
to  tangent  to  rotor  periphery 


A       property  difference 
6       deviation  angle 

e       flow  ratio  =    (Qbe/Qtp) 
p       fluid  density 

<i>       flow  coefficient  =    (C  /U) 
T  m 

2 

\b  work  coefficient  =  (g  Ah  /U  ) 

r  CO 

2 

i|; '       head  coefficient  =  (gAH  /U  ) 


SUBSCRIPTS 

0  total  (static  plus  dynamic)  or  stagnation  value 
of  property 

1  plane  at  entrance  to  rotor 

2  plane  at  exit  from  rotor 

3  plane  at  entrance  to  pump 
be  best-efficiency  point 

f  wall-friction  component 

L  liquid 

m  meridional  or  radial  component  of  velocity 

9  tangential  component  of  velocity 

s  flow-separation  component 

sp  single-phase  component 

T  total  flow,  liquid  plus  vapor 

tp  two-phase  flow 

th  theoretical  or  ideal  conditions 

v  vapor 


INTRODUCTION 

BACKGROUND 

In  a  nuclear  power  plant  of  the  pressurized-water- 
reactor  type  (PWR) ,  sub-cooled  water  at  high  pressure  and 
temperature  is  circulated  through  the  reactor  core  by  several 
main-coolant  pumps,  each  in  its  own  coolant  leg.   Safety 
analyses  of  postulated  loss  of  coolant  accidents  (LOCA) 
require  models  of  the  flow  through  the  main-coolant  pumps 
for  any  one  of  several  operating  modes  in  both  single-  and 
two-phase  flow.   In  an  effort  to  more  precisely  predict  the 
contribution  of  the  main-coolant  pumps  to  flow  through  the 
reactor  core  during  a  LOCA,  there  has  been  a  continuing 
effort  to  improve  these  models,  particularly  for  two-phase 
flow. 

A  semi-empirical  method  for  predicting  the  two-phase 
flow  performance  of  a  centrifugal  pump  was  proposed  by 
Mikielewicz  and  Wilson,   and  was  applied  to  normal  pump 
operation  (forward  rotation  and  flow)  with  promising  results. 
Details  of  their  procedure  are  summarized  below. 

Starting  with  the  Euler  equation  a  theoretically- 
ideal  head-versus-f low  relationship  is  derived  for 
both  single-  and  two-phase  flow. 

By  examining  the  hydraulic  losses  in  two-phase 
flow  it  is  shown  that  for  a  given  pump  the  relative 
pressure  losses — the  difference  between  theoretical 
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and  actual  pressure  rise — are  related  to  the  relative 
pressure  losses  in  single-phase  flow  by  a  function  of 
void  fraction,  a  property  of  the  two-phase  fluid. 

The  ratio  of  relative  pressure  losses  in  two- 
phase  flow  to  those  in  single-phase  flow  is  defined 
as  the  "Head- loss  ratio". 

The  head-loss  ratio  is  then  calculated  by 
combining  the  derived  theoretical  characteristics 
with  experimental  data  for  single-  and  two-phase  flow. 

For  first  quadrant  operation,  good  correlation  is 
shown  between  the  head-loss  ratio  and  void  fraction, 
and  the  result  can  be  used  to  predict  the  two-phase 
characteristics  of  the  pump. 

It  is  shown  that  the  head-loss  ratio  for  a  given 
pump  is  related  to  other  pumps  by  a  function  of  pump 
geometry;  further  experimental  work  is  required  to 
determine  this  relationship. 

PURPOSE 

The  purpose  of  this  report  is  to  determine  if  the 
procedure  is  applicable  to  the  third  quadrant  of  operation 
(reverse-flow  and  rotation) .   The  steps  of  the  procedure,  as 
outlined,  will  be  modified  only  as  required  by  the  different 
flow  geometry,  and  the  terminology  used  will  be  the  same. 


PROCEDURE 

DERIVATION  OF  THEORETICAL  CHARACTERISTICS 

A  centrifugal  pump  subjected  to  conditions  such  that  its 
normal  flow  and  rotation  directions  are  reversed  will  operate 
in  a  manner  similar  to  a  radial-in-f low  turbine.   The  Euler 
equation  can  be  applied  to  determine  the  characteristics  of 
the  pump  in  this  mode  of  operation: 

9c  iho  =  U1C61  '  U2C02  (1) 

where : 

CQ  e  fluid  tangential  velocity 

g_  =  constant  in  Newton's  Law 
3c 

Ah    =    change    in   total   enthalpy  =   h    ,-h   ~ 
o  ^  CJC  ol      o2 

U   =  rotor  peripheral  velocity 
subscript  1  e  rotor  entrance 
subscript  2  =  rotor  outlet 

The  fluid  tangential  velocity  at  rotor  outlet,  CQ2/ 
could  be  significant  because  conditions  of  flow  and  rotor 
speed  would  be  other  than  optimum  and  because,  in  general, 
a  pump  is  not  designed  to  be  an  efficient  turbine.   For 
these  reasons  the  simplifying  assumption  of  zero  outlet 
swirl  will  not  be  made.   This  general  form  of  the  equation 
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is  exact  and  will  yield  the  actual  enthalpy  change  across 
the  pump,  for  any  flow  condition,  if  the  proper  values  for 
the  velocities  are  used. 

In  the  inlet  region  of  the  pump,  Figure  1,  precise 
flow  angles  and  velocities  are  unknown  and  an  assumption 
must  be  made  relating  CQ1  to  the  flow  rate  entering  the 

9  J. 

scroll.   Using  this  assumed  value  of  Cfl.  the  velocity  tri- 
angle at  the  rotor  entrance,  as  shown  in  Figure  2,  can  be 
solved  for  all  velocities  in  terms  of  the  flow  rate. 

At  the  rotor  outlet,  flow  geometry  is  specified  by  the 
design  of  what  is  normally  the  impeller  inlet  eye.   A  typical 
design  with  swept-back  blading  is  shown  in  Figure  3.   Reverse 
flow  and  rotation  will  produce  a  velocity  diagram  as  shown  in 
Figure  4,  which  can  also  be  solved  for  all  velocities  in 
terms  of  the  flow  rate. 

If  the  reverse  flow  through  the  pump  is  idealized  by 
making  simplifying  assumptions,  the  Euler  equation  can  be  used 
to  approximate  the  theoretical  head  and  flow  characteristics 
of  the  pump.   The  following  flow  conditions  are  assumed. 

1.  That  the  fluid  tangential  velocity  at  rotor  entrance, 
C„.,  is  equal  in  magnitude  to  the  absolute  fluid 
velocity,  C->,  entering  the  scroll. 

2.  That  the  relative  flow  angle  at  rotor  outlet,  62, 
is  constant  for  all  flow  rates  and  conditions. 

3.  That  the  deviation  angle  at  rotor  outlet,  62'    is 
zero,  that  is,  8'=B?. 
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PLANE  3 


FLOW,  Q 


PLANE  1 


PLANE  2 


FIGURE  1  -  TYPICAL  CENTRIFUGAL  PUMP  IN  REVERSE-FLOW  AND 
ROTATION 
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U..  -  rotor  peripheral  velocity 
C,  -  fluid  absolute  velocity 
CQ,-  fluid  tangential  velocity 
C  ,-  fluid  radial  velocity 
w,  -  fluid  relative  velocity 
8,  -  fluid  relative  flow  angle 


FIGURE  2  -  VELOCITY  DIAGRAM  AT  IMPELLER  INLET 
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FIGURE  3  -  TYPICAL  IMPELLER-EYE  WITH  REVERSE-FLOW  AND 
ROTATION 
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FIGURE  4  -  VELOCITY  DIAGRAM  AT  IMPELLER  OUTLET 
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With  the  above  assumptions,  equation  1  becomes 

g   Ah  .,  =  U,CA.  -  U0C  n  (2) 

^c    oth     1  91     2  92 

where 

subscript  th  =  theoretical  or  ideal  conditions 

For  the  specific  case  of  loss-less  incompressible  flow, 
the  theoretical  total  enthalpy  change  is  proportional  to  the 
change  in  theoretical  total  head. 

9c  4>\>th  "  9  A[feth 

where 

AH_  =   change  in  total  head  =  H  , -H  ~ 
g    =  gravitational  acceleration 

and  the  expression  for  the  theoretical  total  head  character- 
istics of  the  pump  is: 


9  AHoth  "  U1C91  "  U2C6  2  (3) 
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Theoretical  single-phase-flow  characteristics 

For  single-phase  flow,  the  velocity  triangles  are  solved 
by  using  the  continuity  equation  and  the  assumptions  made  to 
idealize  the  flow. 

In  the  pump  inlet  region,  it  was  assumed  that 


c       =  c 

^01    ^3 


where 

subscript  3  =  plane  at  inlet  to  pump 


From  continuity 


Q  =  C3A3  =  CmlAl 


where 

Q  =   volume  flow  per  unit  time 

A  =   area  normal  to  the  flow 

subscript  m  E  meridional  or  radial  component  of  velocity 

therefore, 


Sl  =  cmiJr  <4) 
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At  the  rotor  outlet,  from  Figure  4, 


C  o 

C_  =  U.  -   m2 


92     2    tan  B2 


Since 


r2 

U0  =  U,  — 

2     1  rx 


where 
r 


2 

—  =   rotor-outlet-to-mlet   radius-ratio 

ri 


and  by  continuity 


Al 

C    =  C    — — 

m2     ml  A_ 


then, 


r2    Al   C  1 
C92  =  Ul  r^  "  jq     tan  32  (5) 


By  substituting  equations  4  and  5,  equation  3  becomes 


9  AHothsp  =  UlCml  SJ  "  Ul  ^  [U1  Pf  "  5^  il^] 


where 

subscript  sp  =    single-phase  flow 
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2 

Dividing  by  U.  and  defining  the  head  coefficient,  ty',    as 


gAH 


.    •  O 


1_  "? 


and  the  flow  coefficient,  <K  ,  as 


A      -     ml 


equation  3  becomes 


r2  2     Al    r2       Al 
rlthsp     r-i      A_   r-i   A~  tan  p.,    1 


This  expression  is  a  linear  function  of  \\j '    versus  §, 
and  is  similar  to  the  expression  derived  for  first-quadrant 
operation. 

Theoretical  two-phase-flow  characteristics 

In  two-phase  separated  flow  it  is  assumed  that  the 
liquid  and  the  vapor  phases  have  separate  mass  flows  and 
that  their  velocity  components  may  also  be  different.   In 
this  case,  equation  2  becomes 


♦ithtp  =  ^[(1-xi)ceLi+xicevi]  -  ;p-[(1-x2)ceL2+x2cev2] 
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where 

x  =  quality  =  m  /mL 
and  m  =  mass  flow  per  unit  time 
subscripts  L  =  liquid 

T  =  total  flow,  liquid  +  vapor 
tpE  two-phase  flow 
v  =   vapor 

From  the  velocity  triangles  at  rotor  inlet  and  outlet 
and  from  the  assumption  that  CQ ,  =  C.,, 

A  A 

^'..    =  i-[(l-X.)C  _  -   —•  +  XnC    .  zr^] 
^lthtp      U]  1    mLl  A"3        1  mVl   A"3 

2 
UZ.  C  C 

2  r  / -i    \  /t   mL2      x  ,    ,,   mv2      «  -, 

"  ^2[(1-x2'  (1-U2  tan  62>+X2(1-U2  tan  B2>  ] 
From  continuity 


mL      mT  (1-x) 
CmL  =  *TK   =  PL(1"a)AT 


m      m_  x 

e   =   v   - 


v    P   A     p   a  Am 
v   v    v    T 


where 

A 
a  =  void  fraction  (•%— ) 

AT 
p  =  fluid  density 
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and  by  defining  the  two-phase-flow  coefficient,  <j>   ,  as 
0     It?  .   _^ 


tp"   A  U     p    A  U 


where 


p    =  mean  density  of  two-phase  flow 
=  (l-a)pL  +apv 


Equation  2  becomes 


♦ithtp  "  jq  ♦tpiI(1-*i)2+«1-*i>2  TT^T  ^  + 


X 


2  (1-al}   ^Ll     2,    ^L  *tP2   n+^2 p_v2, 

1   ax    p     X1J  "  tf^-tan   B2U  (l-a2)  PL2J 


[  (1-x,)  2+— -^-  ^  x2]  }  (7) 

2      a2    pv2   2 


By  defining  a  two-phase-flow  function,  a, 


p 

a  =  (■= )  — 

1-ot   pL 


and  a  two-phase  slip-velocity  ratio,  s, 


C 
v 
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The  quality  of  the  fluid  can  be  expressed  as 


as 
x  = 


1+as 


By  substituting  these  expressions,  equation  7  becomes 

A.  (l+ai)  (1+a^2)  rj2 

*lthtp"   A3    *tpl  (1+ai8i,2  "     (ri} 

<J>      ,         (1+a    )  (l+a0s2) 

[l-tSTfi 2  ]  (8) 

tan    62       (l+a2s2)2 

Applying  the  continuity  equation  to  the  two-phase  flow 
coefficients  at  rotor  inlet  and  outlet  yields, 


ptPi  Ai  ui 

*tp2    "tpl   ptp2   A2   U2 


which  reduces  equation  8  to 


An  r~>  n    *j.  i  P,  i  A,  U, 

... ,      .  _1  ,     f       ,f_2.  2  n   tpl    ,       tpl   1  _1, 

^lthtp  "  A3  9tpl   tpl    {ri'     LX"tan  62   tp2   pfc  2  A2  U2J 


where 


g  d+a)  (1+as2) 
tp     (1+as)2 
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Rearranging: 

,  •         tlh 2   +  r^i  f     +  li  PtPx   ^1   ftp2   1 
vlthtp     ^r]/     LA3  rtpl    rx  Ptp2  A2   tan  62J  ^tpl 

(9) 

If  the  flow  becomes  single-phase,  equation  9  reduces  to 

the  single-phase  case,  equation  6,  because  f   ,,  f   2  an<^ 

p.  ,/p,  2  all  equal  unity  when  the  void  fraction,  a, 

equals  zero. 

As  derived,  equation  9  is  cumbersome  and  requires  a 

knowledge  of  the  relationships  between  inlet  and  outlet 

conditions  to  calculate  values  of  the  head  coefficient  for 

various  conditions  of  void  fraction  and  flow.   By  assuming 

that  fluid  properties  at  rotor  inlet  and  outlet  are  equal, 

equation  9  reduces  to 


r2  2    Al    r?     Al 

i>Wh^   =  "(—)+  [^  +  —  t *.*Z    o    1  f^-,  4*^1      (10) 

lthtp     r-,      A_    r-i  A^  tan  3^    tpl   tpl 


The  assumption  is  valid  in  this  case  because  the  moment  of 

momentum  at  rotor  outlet,  U  C   ,  is  very  small  in  comparison 

2  8  *• 

to  U,C  ,,  and  the  outlet  properties  appear  in  the  term  derived 
from  U2C  2;  thus  the  effect  of  changes  in  fluid  properties 
from  inlet  to  outlet  is  minimal. 
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The  two-phase  flow  function,  f   ,  was  defined  by 

up 


equation  8a  as 


f      (1+a) (1+as2) 


*?"    (1+as)2 


where 


a  =  (J*-)    ^ 

1-a   pL 

C 


For  a  steam-water  mixture,  liquid  and  vapor  densities 
are  given  by  steam  tables  for  various  pressures.   The  slip- 
velocity  ratio,  s,  for  two-phase  flow  in  a  horizontal  pipe, 

2 
was  correlated  as  a  function  of  pressure  and  void  fraction 

as  shown  in  Figure  5.   A  correlation  more  applicable  to  two- 
phase  flow  in  centrifugal  pumps  was  not  available;  therefore, 
the  relationship  described,  though  for  vastly  different 
conditions  of  flow,  was  applied. 

Values  of  the  two-phase  flow  function,  f   ,  can  be 
calculated  as  a  function  of  void  fraction,  a,  for  various 
pressures,  as  shown  in  Figure  6. 

The  theoretical  characteristics  for  third-quadrant 
operation  of  a  centrifugal  pump  in  both  single-  and  two- 
phase  flow  can  be  represented  as  shown  in  Figure  7.   In 
single-phase  flow,  f    is  unity;  in  two-phase  flow,  f    is 

determined  by  pressure  and  void  fraction. 
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10.0 


s=- 


0.5 


0.1 


30  40  '.r.   60  80  100 


200      400  600   1000 
PRESSURE,  psia 


2000 


FIGURE  5  -  SLIP-VELOCITY  RATIO,  s,  VERSUS  PRESSURE 

(from  reference  2) 
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1.20 


1.0 


VOID  FRACTION,  a 


FIGURE  6  -  TWO-PHASE  FLOW  FUNCTION,  f^  ,  VERSUS  VOID 

tp' 

FRACTION,  a,  FOR  200  AND  500  psia 
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♦ith 


ri 


Ai  r?   Ai 

SLOPE  =  -i+  * 


A3  r^  tanB, 


ftPl   ♦! 


FIGURE  7  -  THEORETICAL  HEAD- VERS US -FLOW  CHARACTERISTICS 


27 


FLOW  LOSSES 

At  a  given  rate  of  flow  through  a  centrifugal  pump, 
the  difference  between  the  measured  head  and  the  theoretically 
ideal  head  is  the  result  of  flow  losses  generated  within  the 
pump,  and  changes  in  the  exit  deviation  angle ,  <5  ~  /  from  the 
assumed  ideal  direction.   For  third-quadrant  operation,  in 
single-phase  flow,  this  difference  is  shown  in  Figure  8. 

The  curve  of  actual  head  versus  flow  is  typical  for  this 

3  4 
quadrant,  '   and  the  theoretical  relationship  is  approxi- 
mated from  the  previously  derived  expression  (equation  6) . 

The  actual  head,  AH  ,  is  a  function  of  the  flow-rate,  Q, 

o 

squared;   therefore,  with  flow  increasing  from  zero,  the 
losses  decrease  as  the  square  of  the  flow,  pass  through  a 
minimum,  then  increase  as  the  square  of  the  flow.   The  flow- 
rate  corresponding  to  the  point  of  minimum  losses  is  de- 
fined as  the  best  efficiency  flow,  Q,  ,  for  third-quadrant 
operation. 

In  the  first-quadrant  of  operation,  the  flow  losses 

a 
behave  inAsimilar  manner  and  the  expressions  derived  to 

relate  the  losses  in  two-phase  flow  to  the  losses  in  single- 
phase  flow  may  be  applied  to  the  third-quadrant.   The 
constants  and  flow-rate  ratio  are  those  applicable  to  the 
third  quadrant. 


28 


,  gAHo 


•  \ 


FIGURE  8  -  TYPICAL  THEORETICAL  AND  ACTUAL  PUMP 
CHARACTERISTICS 
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The  flow  losses  are  expressed  as  the  sum  of  wall-friction 
losses  and  flow-separation  losses, 


Ap   =  Ap  .  +  Ap  (11) 

*o    'of    ^os 


where 

Ap   5  loss  of  total,  or  stagnation,  pressure 
subscripts  f  =  friction  component 

s  =   flow-separation  component 

The  wall-friction  component  of  the  losses  in  two-phase 
flow  is  given  by 


Ap  .  .  =  R  K  (-^)2  Q^  (12) 

otpf        pL     tp 


where 


5   1-x  2 
R  =  friction  factor  multiplier   [■= ] 

K  =   constant  at  high  Reynolds  numbers;  includes  pump 

geometry,  the  single-phase  friction  factor,  and 

liquid  density 


By  approximating  the  flow-separation  losses  as  proportional 
to  the  square  of  the  difference  between  the  existing  flow 
and  the  flow  at  best  efficiency  point,   Q,  ,  the  two-phase 
flow-separation  losses  become 
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Ap     m   K,  Q2   [(l-a)(i~-)  -  e]2  (13) 

^otps     4   tp        1+sa  v   ' 


where 


K 


.    =   constant,  depends  on  flow  regime 


e   -  Qbe/QtP 


Adding  these  components  yields  the  expression  for  the 

total  losses  in  two-phase  flow,  Ap    , 

^  *otp 


P 


^H)2  o2    +  k.  o2  ra-awJis.- 


APotP  "  R  K  (p^}   QtP  +  K4  Qtp[(1"a) (ITsT)  "  £]     (14) 


substituting 


tp    ..l-a-,  ,1+a  , 
pT   '   ll-xJ  Ll+saJ 

Li 


The  total  losses  in  two-phase  flow  become 


APotp  =  V'rair'2  +  K6[(1-a»  fe»-£l2}  QtP  (15) 


where 

K-  and  Kfi  are  constants  which  include  geometric  and 
flow-regime  factors 

for  single-phase  flow,  equation  15  reduces  to 


Aposp  ■  V1  +  *6<l-e>2]  QL  <16) 
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The  comparison  between  single-phase  and  two-phase  flow 
is  made  at  like  values  of  flow-rate;  therefore,  the  ratio 
of  two-phase  to  single-phase  losses  is 


APptP_  (iH?>2  *  V(1-a)(£t?>  -  £'2  (17) 

Aposp  1   +  Kc  (1-e)2 

D 


HEAD-LOSS  RATIO 

The  total  losses,  Ap  ,  can  be  expressed  as  the  difference 

o 

between  the  actual  and  the  theoretical  heads  at  a  given  flow- 
rate, 


Ap   =  (AH   -  AH  .,  )  p2-  (18) 

co  o     o  th    g 


In  terms  of  head,  equation  17  becomes 


AH  .   -  AH  ...  pT  (1±|^)2+k  r(l-a)  (^fr-)-e]2 

o  tp o  thtp  _  L     1+sa 6 1+sa n  9) 

AH     -  AH  ..    "  pt  -  ,„  ,,   .  2              K 

osp     othsp  Ktp  1+Kg  (1-e) 


substituting 


L     ,l-x>  ,1+a  . 
p    '   ll-a' ll+sa; 


,     sa 
and  x  =  T+il 
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(20) 


Equation  19  becomes 


AHosp  "  AHothsp     (1-a)  (1+a)  [1+Kgd-e)2] 


For  the  first  quadrant,  the  head-loss  ratio,  H* ,  was 
defined  as 


r*  =    othtp      otp 

"  AH_  .    -  AH 

othsp     osp 


To  prevent  the  introduction  of  an  unnecessary  change, 
the  first  quadrant  definition  will  be  used  in  the  third 
quadrant  as  it  is  numerically  equal  to  the  ratio  in 
equation  20.   Thus,  for  third  quadrant  operation 


othsp     osp 


In  this  functional  relationship:   (1+a/l+sa)  is  a 
function  of  void  fraction  and  pressure;  Kfi  includes  the 
geometry  of  a  given  pump  and  the  flow  regime  of  the  two- 
phase  fluid  (a  function  of  void  fraction) ;  and  e    is  the 
ratio  of  best  efficiency  flow  in  third-quadrant  operation 
to  the  existing  flow.   For  a  given  pump,  operating  at  a 
particular  pressure,  the  ratio  of  two-phase  to  single-phase 
losses  at  like  values  of  e  is  therefore  principally  a 

function  of  the  void  fraction  of  the  two-phase  fluid. 
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Substituting  the  definition  for  the  head  coefficient, 
gAH 
the  head-loss  ratio  can  be  expressed  as 


H*  =  *thtp  "  »tp  (22) 

ythsp    rsp 


CALCULATION  OF  THE  HEAD-LOSS  RATIO 

The  head-loss  ratio  as  a  function  of  void  fraction  for 
third-quadrant  operation  of  a  particular  centrifugal  pump 
can  be  constructed  using  the  derived  theoretical  character- 
istics and  the  actual  characteristics  obtained  from  experi- 
mental data.   The  relationship  between  H*  and  a  for  this 
pump  can  then  be  used  to  predict  the  third-quadrant,  two- 
phase  performance  of  other  pumps  as  single-phase  character- 
istics are  usually  known  or  can  be  easily  determined. 

D.J.  Olson  conducted  experiments  on  a  semi-scale  pump 

7 
of  known  geometry,   in  single-  and  two-phase  flow,  and  in 

the  operating  quadrants  of  interest.   Using  the  results  of 
his  experiments  in  the  third  quadrant,  the  actual  character- 
istics of  the  pump  can  be  calculated.   Information  obtained 

g 

from  the  pump  manufacturer   provided  the  required  constants 

for  calculation  of  the  theoretical  characteristics  for  the 

pump  using  equation  10. 
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The  functional  relationship  between  head-loss  ratio 
and  void  fraction  for  this  particular  pump  is  developed  as 
described  in  the  following  sections.   Details  of  Olson's 
apparatus  and  test  conditions  are  shown  in  Appendix  D. 

Theoretical  head-coefficients 

Details  of  the  semi-scale  pump  used  by  Olson  are  shown 
in  Figures  9  and  10.   With  proper  values  substituted  for  the 
constants,  equation  10  becomes 


*ith  =  -°-1896  +  9-155  ftPi  *i 


As  previously  described,  this  equation  may  be  used  to 
determine  the  theoretical  characteristics  in  both  single- 
and  two-phase  flow. 

Actual  single-phase  head-coefficient 

A  non-dimensional  head-flow  curve  for  the  semi-scale 
pump  can  be  constructed,  as  shown  in  Figure  11,  using  the 
single-phase  data.   Curve-fitting  the  points  by  regression, 
using  the  method  of  least  squares,  this  curve  is  accurately 
approximated  by  the  expression 


iK'    =101.28  (K1,8756  +  0.4 
lsp  1 
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A3    ■    0.0089    ft4 


A^O.069    ft' 
1^=3.875    in 


A2=0.051   ft 
r  =1.6875   in 


FIGURE  9  -  DETAILS  OF  THE  SEMI-SCALE  PUMP 
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FIGURE  10  -  DETAILS  OF  THE  SEMI-SCALE  PUMP  IMPELLER  EYE 
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FIGURE  11  -  HEAD- VERSUS -FLOW  CHARACTERISTICS  OF  SEMI-SCALE 


PUMP  IN  SINGLE-PHASE,  THIRD-QUADRANT  OPERATION 
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Head-loss  ratio 

The  head-loss  ratio  is  calculated  for  each  test 
conducted  in  two-phase  flow,  third-quadrant  operation. 

RESULTS 

The  results  of  these  calculations  are  shown  in  Table  1, 
and  the  relationship  between  head-loss  ratio  and  void 
fraction  for  the  semi-scale  pump  is  shown  in  Figure  12. 
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H*  =   Ithtp    yltp 
^ithsp  *  *iap 
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FIGURE  12  -  HEAD- LOSS  RATIO,  H* ,  VERSUS  INLET  VOID  FRACTION,  a. 
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TABLE  1 
CALCULATION  OF  HEAD-LOSS  RATIO,  H* 


INLET 

VOID 
FRACTION 

ftPi 

*tpl 

*lthtp 

*itP 

^ithsp 

rlsp 

H* 

0.461 

1.022 

0.030 

0.095 

0.919 

0.089 

0.545 

1.81 

0.464 

1.015 

0.018 

-0.022 

0.876 

-0.025 

0.454 

1.88 

0.466 

1.022 

0.010 

-0.094 

0.816 

-0.096 

0.419 

1.77 

0 

1.0 

0.009 

-0.102 

0.469 

-0.102 

0.417 

1.10 

0 

1.0 

0.010 

-0.101 

0.481 

-0.101 

0.417 

1.12 

0 

1.0 

0.024 

0.027 

0.570 

0.027 

0.491 

1.17 

0.241 

1.015 

0.044 

0.215 

0.824 

0.209 

0.683 

1.28 

0.116 

1.009 

0.047 

0.247 

0.811 

0.244 

0.732 

1.16 

0.245 

1.007 

0.044 

0.214 

0.848 

0.211 

0.686 

1.33 

0.056 

1.004 

0.025 

0.043 

0.629 

0.042 

0.502 

1.27 

0.512 

1.067 

0.037 

0.168 

1.004 

0.146 

0.605 

1.82 

0.537 

1.023 

0.022 

0.012 

0.956 

0.008 

0.476 

2.02 

0.631 

1.033 

0.015 

-0.048 

1.184 

-0.052 

0.438 

2.51 

0.401 

1.014 

0.016 

-0.043 

0.788 

-0.046 

0.442 

1.70 
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DISCUSSION 

Figure  12  shows  that  within  the  range  of  variables  of 
Olson's  experiments,  the  head-loss  ratio  of  the  semi-scale 
pump  in  third-quadrant  operation  correlates  very  well  as  a 
function  of  inlet  void  fraction.   An  average  curve  drawn 
through  the  points  can  be  used  to  predict  the  two-phase 
performance  at  other  values  of  void  fraction,  flow,  and 
pressure,  up  to  the  maximum  void  fraction  tested. 

At  low  values  of  void  fraction  (up  to  20%)  the  flow  is 
largely  homogeneous  and  the  losses  in  two-phase  flow  are 
approximately  equal  to  those  in  single-phase.   As  void 
fraction  increases  above  20%,  the  flow  becomes  separated 
and  the  relative  losses  increase  sharply,  until  at  the 
maximum  void  fraction  tested  (6  3%) ,  they  are  more  than 
twice  the  single-phase  losses.   By  the  definition  of  the 
head-loss  ratio,  its  value  is  unity  at  a  void  fraction  of 
100%;  the  characteristics  between  63tand  100%  void-fraction 
are  unknown,  but  it  is  reasonable  to  assume  that  the  head- 
loss  ratio  will  continue  to  increase  then  drop  sharply  to 
unity  at  the  very  high  values  of  void  fraction. 

The  effect  of  pressure  on  the  correlation  of  the  head- 
loss  ratio  is  not  apparent  for  this  limited  amount  of  data. 
Two  series  of  two-phase  tests  were  conducted,  one  at  200  psia 
system  pressure,  the  other  at  500  psia.   The  functional 
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relationship  for  the  head- loss  ratio,  equation  21,  indicated 
that  H*  was  dependent  on  void  fraction,  a,  and  pressure 
(for  a  given  pump) .   In  this  equation,  system  pressure 
affects  the  two-phase  flow  function,  a,  through  its  effect 
on  density  and  the  slip-velocity  ratio,  s.   Both  of  these 
effects  are  overshadowed  by  changes  in  void  fraction,  a. 

OTHER  RESULTS 

The  head-loss  ratio  was  calculated  using  expressions 
other  than  equation  10  for  the  theoretical  two-phase  head 
coefficient. 

The  experimental  data  for  the  fluid  properties  at 
pump-outlet  were  used  to  evaluate  equation  9  and  determine 
the  head-loss  ratio.   Results  are  shown  in  Table  El  and 
Figure  El.   These  results  are  not  significantly  different 
from  the  results  presented  in  Figure  12,  indicating  that 
the  effect  of  the  outlet  tangential  velocity,  Cfi2,  is 
minimal,  as  assumed.   Equation  10  can  be  applied  to  any 
pump,  without  prior  knowledge  of  the  relationships  between 
fluid  properties  at  inlet  and  exit;  therefore,  it  is  much 
more  useful  as  a  means  of  predicting  two-phase  characteristics 

Results  obtained  by  assuming  zero  outlet-swirl  are 

shown  in  Figure  E2,  and  also  show  good  correlation  between 

H*  and  a,  but  indicate  much  higher  losses.   The  procedure  is 

simplified  by  this  assumption,  and  further  experiments  may 

prove  it  to  be  viable  for  certain  pump  designs,  particularly 

those  of  low  specific  speed. 
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APPLICATION  OF  THE  PROCEDURE 

The  two-phase ,  third -quadrant  head-versus-f low 
characteristics  of  a  centrifugal  pump  can  be  predicted,  for 
various  values  of  void  fraction  and  pressure,  using  the  pro- 
cedure outlined  below.   The  following  information  concerning 
the  pump  must  be  known. 

Flow  area  at  the  pump  discharge  flange,  A_. 

Flow  area  of  the  outer  radius  of  the  impeller,  A, . 

Flow  area  at  the  impeller  eye,  A„. 

The  impeller  radius  ratio,  t~/t. . 

Blade  angle  at  impeller  eye,  $1,    as  defined  by 

Figure  3. 
Sufficient  single-phase  data  to  construct  a  plot 
of  head-coefficient  versus  flow-coefficient. 
The  actual  single-phase,  and  the  theoretical  single- 
and  two-phase  head-versus-f low  characteristics  are  determined 
by  constructing  a  graph  as  shown  in  Figure  13.   The  three 

head-coefficients,  <Ksp>  ^ithsp'  and  ^lthtp'  can  then  be 
found  for  particular  values  of  flow-coefficient,  <£,,  and 
two-phase  flow  function,  ftDi* 

At  the  desired  value  of  void  fraction,  a,,  and  system 
pressure,  the  head-loss  ratio,  H*,  is  determined  from 
Figure  12,  and  the  two-phase  flow  function,  f.  -.,  calculated. 
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FIGURE  13  -  ACTUAL  SINGLE-PHASE  AND  THEORETICAL  SINGLE- 


AND  TWO-PHASE  PERFORMANCE 
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The  predicted  values  of  actual  two-phase  head  coeffi- 
cient, ^'   ,  are  calculated  as  a  function  of  flow-coefficient, 
(J),  ,  from 


rltp    ylthtp      irlthsp   ylsp 


Experiments  on  pumps  of  different  geometries  will  be 
required  to  determine  the  effect  of  geometrical  factors  such 
as  specific  speed,  N  ,  on  the  functional  relationship  of  the 
head-loss  ratio,  H*  (equation  21) .   Until  such  data  becomes 
available,  this  procedure  may  be  used  to  determine  approximate 
two-phase  characteristics. 
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APPENDIX  A 

DETAILS  OF  THE  DERIVATION  OF  THE  THEORETICAL 

TWO-PHASE  HEAD  COEFFICIENT 


In  two-phase  separated  flow,  the  theoretical  two-phase 
head  coefficient  is  given  by 


♦ithtp  =  ff7!(1-Xl)C9Ll+xlCevl1-r|[(1-x2'C6L2+x2C9v2) 

1  (Al) 

In  general,  velocities  of  the  two  phases  are  given  by 


m_      m_(l-x) 


L    PL  AL    PL(1"a)AT 


m       m_  x 
c   =   v        T 


v    p   A     pa  Am 
^v   v    Kv     T 


Assuming  C,,  =  C_ 


m  (1-x.) 

=  C, 


LI     L3    pL1(l-a1)A3 


mm  x.. 

=  c     T  1 


vl     v3    Pvl  a1  A3 


By  continuity 


C  T1  pT  ,  (1-cu)  A    C   ,  p  ,  an  A 
mLl   LI      1    1  _  mvl   vl   l   1 
mm  =  constant  =  ? - 
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Then 


c    _  Cnil  PU(1"al)Al   (1-X1} 


6L1        (l-x^         pL1(l-a1)A3 


and 

Al 
c    =  c    — =■ 
8L1     mLl  A 

Similarly 


Al 
C     =  C     — — 
0vl    mvl  A 


Substituting  into  equation  Al 


mLl    1  ,      mvl   1 


*ithtP  =  (1-xi'—  a;  +  xi  —  A 


[(1-'2»;2ceL2  +  Wew, 

ul  ul 


From  the  general  expressions  for  velocity 


mT  (1-x) 
C. 


mLl    pL1(l-a1)A1 


and 


mT  Xl 


'mvl    pvl  a1  A1 
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Substituting  these  expressions,  the  first  two  terms 
on  the  right  hand  side  of  equation  A2  become 

2  •  2  • 

(l-x1)   mT    Ax    x±   mT       Ax 


U1PL1(l-a1)A1  A3    pvl  al  U1A1  A3 


Rearranging  and  multiplying  by 


(l-o^)  pL1  +  a1  pvl 
(1-c^)  pLl  +  ax  pvl 


where 


(l-a)pL  +  apv  e  ptp 


They  become 


2         2 

m  A  (1-x^        x 

Substituting  the  definition  of  the  two-phase  flow  coefficient 


,1    -  ^ 


fcPv    ptPl  Al  Ul 


the  two  terms  become 


Al  A     m     2    (!-«!> 2  «i  Pvi   *1  PL1  (1-al>  ..  ..2. 
A7*tpl  '"-i)  +^7TT^J pv1  a, 
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by  definition,  the  two-phase  flow  function,  a,  is 


a  =  [JL.]  !z 

1-a   pL 


and  the  quality,  x,  becomes 


as 
x  = 


1+as 


where 


C 
v 


substituting  these  expressions,  the  first  two  terms  become 


2 
A,  2aisi      (ais-|)  i 

jt  ^toi  C1+ai]  [1-i+i-i-  +  — — 2  (1+  ? ] 

A3   tpl     X     l+^i3!    l+2a1s1+(a1s1)2     a 
performing  the  indicated  operations  and  rearranging  yields 


A±  (1+a^  (l+a1s2) 

<J>. 


A3  ^    (1+a^,2 


Substituting  this  expression,  equation  A2  becomes 


Al       .       r(1-X2)U2C9L2  +  X2U2C8v21 

*lthtp  ~  aJ  ^tpi  rtPl        u2        u2 
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where 


(l+a1)  (1+a.^) 


tpl  "    z,^     ,2 
^       (l+a1s1) 

From  the  outlet  velocity  triangle.  Figure  4  ,  the  last  two 
terms  in  equation  A2  reduce  to 


_  ftp2  ({)tp21  52,2 
l±    tan  62   J  ^r 


using  the  same  procedure  applied  to  the  first  two  terms, 
Equation  A2  becomes 


A!  a     4=      ,^2  n    ftp2  4'tp21 
*lthtp  "  AT  ^tpl  ftpl  -(rY}   [1  "   tan  B2   ] 


The  two-phase  flow  coefficients,  4>tD'  at  inlet  and  outlet 
can  be  related  by  the  continuity  equation, 


A      -  *      ptPi  Ai  °i 

*tp2    *tpl  ptp2   A2   U2 


Substituting  this  expression  and  rearranging  yields 

r~  -,    A,         r_  P.^t  A   f.^0 
^lthtp  '  "ir1;     LA3  rtpl   r   ptp2  A2  tan  82J  vtpl 

(A3) 
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APPENDIX  B 
DETAILS  OF  THE  CALCULATION  OF  HEAD-LOSS  RATIO,  H* 

The  head-loss  ratio,  H*,  for  the  various  values  of  void 
fraction,  a,  was  calculated  using  the  relationships  previously- 
presented.   The  data  from  the  experiments  conducted  by  Olson 
are  included  in  Appendix  C.   Values  for  fluid  densities 
were  interpolated  from  reference  9 . 

Table  Bl  lists  the  single-phase  characteristics  used 
to  determine  the  single-phase  head  coefficient,  tyl       ,  as  a 

JL  O  S"' 

function  of  flow  coefficient,  <J>.  .   Tables  B2  through  B5 
present  the  calculation  of  the  head-loss  ratio,  H* . 

The  test  number  column  refers  to  the  experimental  data. 
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TABLE  Bl 

CALCULATION  OF  THE  SINGLE-PHASE  HEAD  COEFFICIENT,  ty 1 
lsp 

AS  A  FUNCTION  OF  FLOW  COEFFICIENT,  <J> . 


SINGLE- 

gAH 

PHASE 
TEST 

(ft2/s 

NUMBER 

78 

1330 

82 

1459 

86 

164.22 

87 

570 

88 

763 

89 

1346 

92 

74 

83 

4176 

75 

605 

76 

960 

77 

1513 

91 

299.5 

95 

2380 

98 

1204 

150 

4405 

106 

3130 

107 

3916 

108 

4634 

149 

7065 

147 

7467 

148 

9792 

152 

9982 

u. 


ml 


lsP   2 


<J>. 


(ftVsec2)   (ft/sec)   (ft/sec)    (gAH/U")    ^ml/11!5 


13.2 

14.2 

13.4 

13.5 

13.5 

13.6 

13.6 

20.8 

27.1 

27.2 

27.1 

27.04 

53.78 

54.15 

54.25 

74.6 

74.7 

74.6 

74.5 

90.9 

91.2 

92.0 


3.188 

3.275 

0.75 

1.962 

2.308 

3.188 

0 
3.635 
1.255 
2.121 
3.231 

0 

2.568 

0 
4.95 
1.414 
2.640 
3.534 
5.756 
5.02 
6.492 
6.65 


7.64 

7.22 

0.91 

3.15 

4.22 

7.40 

0.4 

9.63 

0.82 

1.30 

2.06 

0.41 

0.82 

0.41 
1.50 
0.56 
0.70 
0.83 
1.27 
0.90 
1.18 
1.18 


0.24 
0.23 
0.06 
0.15 
0.17 
0.23 

0 
0.17 
0.05 
0.08 
0.12 

0 

0.05 

0 
0.09 
0.02 
0.04 
0.05 
0.08 
0.055 
0.07 
0.07 
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TABLE 

B2 

CALCULATION 

OF  TWO-PHASE  FLOW  FUNCTION,  f .   , 

tp 

i 

AND  TWO-PHASE  DENSITY 

,    P.  /  AT 
tp 

PUMP  INLET 

TWO- 

VAPOR 

LIQUID 

VOID 

PHASE 
TEST 
NUMBER 

DENSITY 
\1    3 
lbn/ft 

DENSITY 

PL1  3 
lbn/ft 

FRACTION 
al 

ftPi 

Ptpl 

lb  /ft3 

m 

149 

0.4377 

54.38 

0.548 

1.0299 

24.82 

150 

0.4443 

54.33 

0.461 

1.0218 

29.49 

151 

0.4508 

54.28 

0.464 

1.0151 

29.30 

152 

0.4354 

54.40 

0.466 

1.0218 

29.25 

153 

0.4354 

54.40 

0** 

1.0 

54.40 

154 

0.4420 

54.35 

0** 

1.0 

54.35 

155 

0.4309 

54.43 

0** 

1.0 

54.43 

156 

0.4288 

54.45 

0.241 

1.0150 

41.43 

157 

1.1002 

50.53 

0.116 

1.0086 

44.80 

158 

1.1066 

50.51 

0.245 

1.0069 

38.41 

159 

1.1042 

50.52 

0.056 

1.0039 

47.75 

160 

1.0932 

50.57 

0.512 

1.0671 

25.24 

161 

1.1156 

50.46 

0.537 

1.0232 

23.96 

162 

1.1156 

50.46 

0.631 

1.0327 

19.32 

163 

1.1453 

50.33 

0.401 

1.0143 

30.61 

164 

1.1202 

50.44 

0.647 

1.0348 

18.53 

*Tests  149  through  156  were  conducted  at  nominal  system 
pressure  of  200  psia,  from  Figure  5,  slip  velocity  ratio, 
s  =  2.8.   Tests  157  through  164  were  conducted  at  nominal 
system  pressure  of  500  psia,  s  =  2. 

**Negative  void  fractions  in  data  were  taken  to  be  zero. 
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TABLE  B3 


CALCULATION  OF  TWO-PHASE  FLOW  FUNCTION,  f.  , 

tp 

AND  TWO-PHASE  DENSITY,  pt  ,  AT  PUMP  OUTLET 
tp  


TWO- 

VAPOR 

LIQUID 

VOID 

PHASE 

TEST 

NUMBER 

DENSITY 

V2  3 
lb  /ftJ 
irr 

DENSITY 

PL2  - 

lb  /ftJ 
m 

FRACTION 
a2 

tp2 

PtP2 

lb  /ft 
m 

149 

0.4261 

54.47 

0.840 

1.1070 

9.07 

150 

0.4113 

54.60 

0.867 

1.1234 

7.62 

151 

0.4197 

54.52 

0.781 

1.0767 

12.27 

152 

0.4197 

54.52 

0.558 

1.0298 

24.33 

153 

0.4009 

54.68 

0.008 

1.0002 

54.25 

154 

0.3967 

54.71 

0.007 

1.0002 

54.33 

155 

0.3946 

54.73 

0.730 

1.0568 

15.07 

156 

0.3864 

54.80 

0.789 

1.0742 

11.87 

157 

1.0514 

50.76 

0.702 

1.0405 

15.86 

158 

1.0646 

50.70 

0.756 

1.0509 

13.18 

159 

1.0624 

50.71 

0.428 

1.0147 

29.46 

160 

1.0646 

50.70 

0.696 

1.0400 

16.15 

161 

1.0844 

50.61 

0.912 

1.1065 

5.44 

162 

1.0844 

50.61 

0.988 

1.0860 

1.68 

163 

1.1173 

50.46 

0.744 

1.0505 

13.75 

164 

1.0998 

50.54 

0.834 

1.0736 

9.31 

• 

(1+a)  (1+as2) 

ftp 

(1+as)2 

a   -  (T 

p 
-a   PL 

ptP  =  (1 

-a)  PL  +  ap 

V 
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TABLE  B4 


CALCULATION  OF  TWO-PHASE  FLOW  COEFFICIENT,  <f>.  , 
Ttpl 

AND  TWO-PHASE  HEAD  COEFFICIENT,  \p  ' 


TWO-     FLOW         ROTOR        <J> .  . 

tpl 


PHASE    RATE,  Q      SPEED,  N 

TEST 

NUMBER 


TEST     ,   .  ,    .     .     ,    /  •  x     (0.951  Q/N) 
(gal/mm)    (rev/mm)  ' 


149 

-85.9 

150 

-51.5 

151 

-30.3 

152 

-17.3 

153 

-16.2 

154 

-16.4 

155 

-39.8 

156 

-73.4 

157 

-80.0 

158 

-73.8 

159 

-42.6 

160 

-61.4 

161 

-36.2 

162 

-25.3 

163 

-26.5 

164 

-112.9 

0 

-1608  0.030 

-1600  0.018 

-1607  0.010 

-1604  0.009 

-1602  0.010 

-1600  0.024 

-1603  0.044 

-1607  0.047 

-1603  0.044 

-1600  0.025 

-1594  0.037 

-1597  0.022 

-1600  0.015 

-1603  0.016 
0 


HEAD 

AH 

*itP 

o 

(ft) 

(28,20 

AH 

°) 

2  ' 
N 

40.2 

— 

84.3 

0.919 

79.6 

0.876 

74.8 

0.816 

42.8 

0.469 

43.8 

0.481 

51.8 

0.570 

75.1 

0.824 

74.3 

0.811 

77.3 

0.848 

57.1 

0.629 

90.5 

1.004 

86.5 

0.956 

107.5 

1.184 

71.8 

0.788 

76.4 

— 
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TABLE  B5 

CALCULATION 

OF  HEAD- 

•LOSS  RATIO,  H* 

TWO- 
PHASE 
TEST 
NUMBER 

*ithtP 

vltp 

^lthsp 

'  lsp 

H* 

INLET 

VOID 

FRACTION 

°1 

149 
150 

0.095 

0.919 

0.089 

0.545 

1.81 

0.461 

151 

-0.022 

0.876 

-0.025 

0.454 

1.88 

0.464 

152 

-0.094 

0.816 

-0.096 

0.419 

1.77 

0.466 

153 

-0.102 

0.469 

-0.102 

0.417 

1.10 

0 

154 

-0.101 

0.481 

-0.101 

0.417 

1.12 

0 

155 

0.027 

0.570 

0.027 

0.491 

1.17 

0 

156 

0.215 

C.824 

0.209 

0.683 

1.28 

0.241 

157 

0.247 

0.811 

0.244 

0.732 

1.16 

0.116 

158 

0.214 

0.848 

0.211 

0.686 

1.33 

0.245 

159 

0.043 

0.629 

0.042 

0.502 

1.27 

0.056 

160 

0.168 

1.004 

0.146 

0.605 

1.82 

0.512 

161 

0.012 

0.956 

0.008 

0.476 

2.02 

0.537 

162 

-0.048 

1.184 

-0.052 

0.438 

2.51 

0.631 

163 

-0.043 

0.788 

-0.046 

0.442 

1.70 

0.401 

164 

— 

— 

— 

— 

— 

— 

nth  ■  -°-1896  +  9-155  ftPi  *i 


ibl         =    101.28    <^*8756    +    0.4 
ylsp  Yl 


H*  _  *ithtp  "  HtP 
*ithsp  "  *isP 
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APPENDIX  C 
EXPERIMENTAL  DATA 

Tables  CI  and  C2  are  reproduced  from  reference  7. 

Explanatory  notes  for  single-phase  tests: 

a.  Based  on  turbine  flowmeter  measurements. 

b.  Positive  pump  AP  represents  higher  pressure 
at  discharge  nozzle. 

c.  Based  on  fluid  density  measurement  at  pump 
inlet. 

d.  Fluid  density  at  pump  inlet  (suction  for 
normal  operation,  discharge  for  reverse  flow) 

e.  Inlet  and  outlet  refer  to  locations  with 
respect  to  direction  of  flow. 

f.  Flow  rate  obtained  from  venturi  flowmeter 
in  test  loop. 

Explanatory  notes  for  two-phase  tests: 

a.  Inlet  flow  based  on  loop  mass  flow  rate  and 
test  section  inlet  (suction  or  discharge) 
fluid  density. 

b.  Positive  AP  indicates  higher  pressure  at 
pump  discharge. 

c.  Based  on  inlet  conditions. 

d.  Inlet  and  outlet  refer  to  the  direction  of 
flow. 
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e.  Calculated  using  measured  density  and 
assumed  saturation  conditions. 

f.  Measured  directly. 
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APPENDIX  D 
DESCRIPTION  OF  EXPERIMENTS 

The  semi-scale  pump  used  in  the  experiments  had  the 
following  characteristics: 

Rated  flow  -  180  gpm 
Rated  head  -  192  ft. 
Specific  speed  -  926 
Rated  torque  -  417.6  in-lb 
Details  of  the  pump  are  shown  in  Figures  9  and  10. 
The  test-loop  used  for  single-phase  tests  is  shown  in 
Figure  Dl.   The  two-phase  test-loop  is  shown  in  Figure  D2. 
Superheated  steam  and/or  saturated  water  were  supplied  by 
an  oil-fired  boiler. 

Flow  conditions  imposed  for  third  quadrant  tests  are 
as  follows. 

Single-phase  tests 

Impeller  speed  -  up  to  minus  2700  rpm 
Flow  rate  -  0  to  minus  207  gpm 
Pressure  at  pump  inlet  -  60.7  to  216.6  psig 
Pump  head  -  5.1  to  310.0  ft  of  fluid 

Two-phase  tests 

Impeller  speed  -  0  and  minus  1600  rpm 

Flow  rate  -  minus  16.2  to  minus  112.9  gpm 

System  pressure  -  200  and  500  psia 

Pump  head  -  40.2  to  107.5  ft  of  fluid 
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T      Thermocouple 

P       Pressure  Top 

Differential    Pressure   Transducer 
/l\      Blanking  Plates  Installed  for  Forward  Flow 
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FIGURE  Dl  -  SINGLE-PHASE  TEST  LOOP  (from  reference  7) 
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FIGURE   D2    -   TWO-PHASE   TEST   LOOP    (from  reference    7) 
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APPENDIX  E 
OTHER  RESULTS 

This  section  summarizes,  for  purposes  of  comparison, 
the  correlation  of  H*  vs  a  obtained  using  alternate  methods 
to  compute  the  theoretical  characteristics. 

The  experimental  data  were  used  to  evaluate  equation  9, 
thereby  making  it  semi-empirical  rather  than  theoretical. 
The  head  coefficients  obtained  in  this  manner  were  used  to 
calculate  the  head-loss  ratio.   Results  are  shown  in  Table  El 
and  Figure  El. 

ZERO  SWIRL  AT  OUTLET 

The  theoretical  characteristics  can  be  derived  assuming 
that  outlet  swirl  is  negligible.   In  this  case  the  expression 
for  the  theoretical  head  coefficient,  ty]thf    is 

Al 

*ith  =  a;    ftPi  *i 

Results  are  shown  in  Table  E2,  and  the  correlation  of 
H*  versus  ^  is  shown  in  Figure  E2. 


67 


3.0 

/ 

/ 

o 

2.0 

< 

£ 

H* 

0 

\ 

■ 

1.0 

0 

^200   psii 
O500   psia 

0.2        0.4        0.6 

INLET  VOID  FRACTION,  a. 


O08 


1.0 


FIGURE  El  -  HEAD-LOSS  RATIO  VERSUS  INLET  VOID-FRACTION 


USING  EQUATION  9  FOR  ^thsp 
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TABLE  El 
CALCULATION  OF  THE  HEAD-LOSS  RATIO,  H*  ,  USING 
EXPERIMENTAL  DATA  TO  EVALUATE  EQUATION  9 


TWO- 
PHASE 
TEST 
NUMBER 

^lWse 

H* 

al 

150 

0.199 

1.581 

0.461 

151 

0.006 

1.816 

0.464 

152 

-0.092 

1.763 

0.466 

153 

-0.102 

1.100 

0 

154 

-0.101 

1.124 

0 

155 

0.096 

1.023 

0 

156 

0.336 

1.030 

0.241 

157 

0.341 

0.964 

0.116 

158 

0.306 

1.142 

0.245 

159 

0.059 

1.238 

0.056 

160 

0.189 

1.776 

0.512 

161 

0.097 

1.835 

0.537 

162 

0.128 

2.156 

0.631 

163 

-0.022 

1.660 

0.401 

f'UJ   =  -0.1896  +  (8.1342  f.  ,  +  1.021  -^   f.  9)  <J>.  - 
ylthtp  tpl        ptp2    p     p 


Values  for  ftpl*  ftp2'  ptpl'  and  Ptp2  were  obtained  from 
Appendix  B. 


+ 
where  subscript  se  =  semi-empirical 
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FIGURE  E2  -  HEAD-LOSS  RATIO  VERSUS  INLET  VOID  FRACTION  FOR 


ZERO-SWIRL  AT  OUTLET,  CQ2  =  ° 
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TABLE  E2 
CALCULATION  OF  THE  HEAD-LOSS  RATIO,  H*, 
ASSUMING  ZERO  SWIRL  AT  ROTOR  OUTLET,  CQ2    =  0 


TWO- 
PHASE 
TEST 
NUMBER 

150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 


^ispth 

*ltpth 

H* 

0.248 

0.253 

2.24 

0.146 

0.148 

2.36 

0.083 

0.085 

2.18 

0.078 

0.078 

1.15 

0.079 

0.079 

1.19 

0.193 

0.193 

1.27 

0.354 

0.359 

1.41 

0.385 

0.388 

1.22 

0.356 

0.358 

1.48 

0.206 

0.207 

1.43 

0.298 

0.318 

2.23 

0.175 

0.179 

2.58 

0.122 

0.126 

3.35 

0.128 

0.130 

2.10 
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